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ABSTRACT 

We examine effects on primordial nucleosynthesis from a truly random spatial dis- 
tribution in the baryon-to-photon ratio (77). We generate stochastic fluctuation spectra 
characterized by different spectral indices and root-mean-square fluctuation amplitudes. 
For the first time we explicitly calculate the effects of baryon diffusion on the nucleosyn- 
thesis yields of such stochastic fluctuations. We also consider the collapse instability of 
large-mass-scale inhomogeneities. Our results are generally applicable to any primordial 
mechanism producing fluctuations in rj which can be characterized by a spectral index. 
In particular, these results apply to primordial isocurvature baryon fluctuation (PIB) 
models. The amplitudes of scale-invariant baryon fluctuations are found to be severely 
constrained by primordial nucleosynthesis. However, when the 77 distribution is charac- 
terized by decreasing fluctuation amplitudes with increasing length scale, surprisingly 
large fluctuation amplitudes on the baryon diffusion scale are allowed. 

Subject headings: cosmology: theory - early universe - nuclear reactions, nucleosynthesis, 
abundances - large-scale structure of universe - dark matter 
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1. Introduction 

In this paper we investigate the effects of baryon 
diffusion on the big bang nucleosynthesis abundance 
yields for models which are characterized by a ran- 
domly distributed baryon-to-entropy ratio. Such mod- 
els are distinctively different from the well-studied 
models (Malaney & Mathews 1993) of adiabatic in- 
homogeneities or baryon inhomogeneitics initially in- 
spired by a first-order QCD phase transition. We dis- 
cuss models both with and without the collapse of 
a significant fraction of the large-mass scale, high- 
density regions. Our study is motivated by the 
primordial, isocurvature, baryon (PIB) fluctuation 
model of structure formation (Peebles 1987ab). How- 
ever, the results and constraints derived here apply to 
any theory which produces stochastic fluctuations in 
the baryon-to-entropy ratio which can be character- 
ized by a spectral index. 

Primordial, isocurvature, baryon fluctuations have 
been proposed (Peebles 1987ab) as a model for the 
formation of cosmic large-scale structure. In this 
model it is envisioned that at early times the universe 
was filled with fluctuations in the baryon-to-entropy 
ratio characterized by increasing fluctuation ampli- 
tudes with decreasing mass scale. PIB models with 
primordial power spectra of the right normalization 
and spectral index ( — 1 < n ^ 0, where n — cor- 
responds to white noise) may be quite successful in 
reproducing the statistical properties of the observed 
large-scale structure at the present epoch (Sugino- 
hara & Suto 1992; Cen, Ostriker, & Peebles 1993). 
These models have even certain advantages over other 
popular structure formation models (e.g. cold dark 
matter or mixed dark matter with n w 1) including 
the early formation of structure and the small pecu- 
liar velocity flows on small scales. An often invoked 
feature of PIB models is the collapse of Jeans-mass 
size (Mj - 3 x lO 5 M (f2 6 /i 2 /O.l)), non-linear en- 
tropy fluctuations around the epoch of recombination 
and the early formation of baryonic, collisionless dark 
matter in the form of black holes or brown dwarfs 
(Hogan 1993, Loeb 1993). 

PIB models, however, are severely constrained by 
current observations of the anisotropics in the cos- 
mic microwave background radiation. In contrast to 
cold dark matter or mixed dark matter models, PIB 
models require an early reionization (Peebles 1987a; 
Gorski & Silk 1989; Chiba, Sugiyama, & Suto 1993; 
Hu & Sugiyama 1994). Moreover, this reionization 



should only be partial and occur in the right redshift 
range. It remains to be shown if this can be accom- 
plished self-consistently in such models. 

PIB models usually assume that the fractional con- 
tribution of baryons to the critical density is ~ 
0.1—0.2 (including remnants, as well as diffuse baryons). 
This value is much larger than that allowed by the 
light element abundance constraints of a homoge- 
neous big bang nucleosynthesis, i.e. 0.007 < fi^/i 2 < 0.023, 
where h is the Hubble parameter in units of 100 km 
s _1 Mpc -1 (Smith, Kawano, & Malaney 1993; Copi, 
Schramm, & Turner 1995; Cardall & Fuller 1996). 
This fact, and the fact that the existence of entropy 
fluctuations in general significantly change the abun- 
dance yields during primordial nucleosynthesis, has 
led to at least three independent derivations of nu- 
cleosynthesis constraints on the PIB fluctuations for 
small mass scales (Gnedin, Ostriker, & Rees 1995; 
Jcdamzik & Fuller 1995; Copi, Olive, & Schramm 
1995). On the other hand, it had already been pointed 
out in earlier more schematic studies that the up- 
per limit on f2& may be substantially larger that that 
derived from homogeneous big bang nucleosynthe- 
sis when a significant fraction of high-density regions 
collapse and "swallow" their high 4 He and 7 Li yields 
(Harrison 1968; Zcldovich 1975; Rees 1984; Sale & 
Mathews 1986). 

Probably the most important result of all three re- 
cent investigations of big bang nucleosynthesis with 
stochastic entropy fluctuations and collapse is that 
there is no agreement between calculated primordial 
abundances and observationally inferred abundance 
limits for any Clf, when significant fluctuation power 
is present on mass scales below the Jeans mass. This 
is because high-density fluctuations on scales smaller 
than the Jeans mass don't collapse but rather expand 
and mix their high 4 He and 7 Li yields with the re- 
maining diffuse baryons (Jedamzik & Fuller 1994). 
In other words, a simple extrapolation of the pre- 
ferred stochastic PIB fluctuation spectrum to small 
mass scales with a spectral index (— 1 < n < 0) seems 
to be ruled out by nucleosynthesis considerations. 

In an attempt to circumvent this constraint, Gnedin 
et al. 1995 focussed on the possibility that entropy 
fluctuations are strongly phase correlated. They con- 
sidered super- Jeans-mass size fluctuations described 
individually by a radial power law- or top hat- dis- 
tribution. Retaining the stochastic nature of fluc- 
tuations Jedamzik & Fuller 1995 (hereafter; JF95) 
speculated on the possibility of a fluctuation cut- 
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off on mass scales above the Jeans mass. In both 
cases the authors concluded that the total Of, in rem- 
nants and diffuse baryons may be as large as that 
reflected by dynamical estimates of the density pa- 
rameter, il ~ 0.1 — 0.2 (Peebles 1986; Kaiser et 
al. 1991; Dekel et al. 1992; Strauss et al. 1992), 
without being in conflict with nucleosynthesis abun- 
dance constraints. This is true as long as the high- 
density regions collapse efficiently JF95 also pointed 
out that PIB-like models can produce intrinsic, large- 
mass scale variations in the primordial abundances 
which may possibly be inferred from future observa- 
tions of ( 2 H/H)-ratios in high-rcdshift Lyman-limit 
clouds (see also Jedamzik & Fuller 1996). 

In any case, all previous studies did not adequately 
address the effects of dissipative processes during the 
nucleosynthesis era on the primordial nucleosynthe- 
sis yields. It is well known that baryon diffusion 
and hydrodynamic expansion during primordial nu- 
cleosynthesis affects the abundance yields when en- 
tropy fluctuations on mass scales, M, in the range 
1O~ 22 M < M <10" 10 M Q are present. It is, perhaps, 
a bold extrapolation to infer the amplitudes of en- 
tropy fluctuations on these small mass scales from the 
amplitudes of entropy fluctuations on mass scales rel- 
evant to structure formation (Al ^> Mj). Neverthe- 
less, one would expect that inflation generated spectra 
should extend over many orders of magnitude. And 
in any event, disregarding large-scale structure forma- 
tion, it is instructive to know which kind of stochastic 
fluctuation spectra (existing with appreciable ampli- 
tudes on either only small scales or both, small and 
large scales) are compatible with nucleosynthesis con- 
siderations. 

Gnedin et al. (1995) have attempted to account 
for dissipation during nucleosynthesis by using the 
small-scale, inhomogeneous nucleosynthesis yields de- 
rived by Mathews et al. (1990). The calculations 
by Mathews et al., however, correspond to the abun- 
dance yields produced in a regular lattice of fluctua- 
tion sites, such as might be produced during a first 
order QCD transition. In contrast, the few existing 
baryogenesis scenarios capable of producing PIB-like 
fluctuations suggest a stochastic nature of fluctua- 
tions (Yokoyama & Suto 1991; Dolgov & Silk 1992). 
Moreover, Gnedin et al. used only those fluctua- 
tion parameters which optimize the upper limit on 
fib, which is clearly an unrealistic oversimplification. 
JF95 attempted to qualitatively estimate the effects of 
dissipation during nucleosynthesis on the abundance 



yields for a random entropy distribution. They sug- 
gested that the 4 He and 7 Li yields of a stochastic 
distribution are generally lowered by the effects of 
baryon diffusion but should still be higher than the 
corresponding yields in a homogeneous big bang nu- 
cleosynthesis scenario at the same rib- 
In this paper we compute primordial nucleosyn- 
thesis yields for random entropy fluctuations with 
different spectral indices and fluctuation amplitudes. 
We explicitly incorporate effects of baryon diffusion 
during the nucleosynthesis epoch. Our calculations 
serve to place constraints on the homogeneity in the 
baryon-to-photon ratio independent of any gravita- 
tional collapse. We also combine our results with PIB- 
like models of structure formation in order to derive 
approximate constraints on the allowed isocurvature 
power spectrum. In Section 2 we give the relevant 
physical scales entering the problem, introduce our 
models of stochastic baryon inhomogeneity which are 
based on the models of JF95, and explain our nu- 
merical techniques for the calculation of primordial 
abundance yields. In Section 3 we present our results 
and in Section 4 we draw conclusions. 

2. Random Baryon Fluctuations and Primor- 
dial Nucleosynthesis 

2.1. Physical Length and Mass Scales 

In what follows we will give all length scales in 
comoving units at the present epoch. Often we will 
associate a length scale, A , with a baryonic mass 
scale, Mh(A), by employing 

M b = Jp b A 3 , (1) 
6 

where p~b is the average baryon density at the present 
epoch. There are three relevant physical scales in our 
problem. First, there is the neutron diffusion length 
scale, d n (T w 500keV) , determined at the epoch 
of weak freeze-out in the early universe. This scale 
divides length scales, L , which have been homoge- 
nized by baryon diffusion before the onset of primor- 
dial nucleosynthesis, L < d n , from length scales on 
which inhomogeneities in the baryon-to-photon ratio 
r\ can persist through the beginning of primordial nu- 
cleosynthesis, L > d n , 

d„(T« 500keV)«4x 10~ 5 pc . (2) 

In deriving Eq. || we assume flbh 2 < 1 locally (see 
Jedamzik & Fuller 1994). Using Eq. [| with d n = A, 
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we find the corresponding mass scale 

Msookcv « 1 x 1CT 22 M Q (^ 2 /0.0125) . (3) 

The second length scale is the neutron diffusion 
scale at the approximate completion of freeze-out 
from nuclear statistical equilibrium 

d n (T w lOkeV) » 8 x 10" 2 pc(fi b /i 2 /0.0125)"5 . (4) 

This scale corresponds to a baryonic mass scale 

^fokcv « 1 x lO" 12 M (fi fc /i 2 /O.O125)-* . (5) 

We note here that the relevant M d at the comple- 
tion of the freeze-out process is also sensitive to the 
local flbh 2 because the termination temperature of 
the primordial nucleosynthesis process is itself a func- 
tion of fifc/i 2 . Small local flbh 2 imply relatively large 
M d 3> 1O _12 M . Regions separated by length scales, 
L > rfiokcV are not m con tact through baryon dif- 
fusion during the epoch of primordial nucleosynthe- 
sis. If any one region of the universe is homogeneous 
on length scales ^xokeV but inhomogeneous on larger 
scales, the nucleosynthesis yields can be simply com- 
puted by an appropriate abundance average over dif- 
ferent homogeneous big bang nucleosynthesis scenar- 
ios at different rj. This is what has been done in all 
prior studies (Gnedin et al. 1995; JF95; Copi et 
al. 1995). However, when small-scale inhomogeneities 
are present a reliable calculation has to at least cover 
the length scale range from d5 00koV to d" oke y Our 
calculations achieve this. Finally, there is the Jeans 
length 

A. 7 «5xlO- 2 Mpc(^) 2 , (6) 
yielding the Jeans mass 

M b j » 3 x lO 5 M (fi b /i 2 /O.l)"^ . (7) 

The ultimate fate of overdense regions with M > M h j 
is collapse, whereas isolated, overdense regions with 
M < Mj will expand and mix their nucleosynthesis 
yields with the remaining diffuse baryons. It should 
be noted that Mj is not the true Jeans mass before the 
epoch of recombination. It rather compares the effects 
of baryonic pressure with baryonic self-gravity and 
determines if a fluctuation will expand or contract. 



It coincides with the post-recombination Jeans mass 
when all photon stresses can be neglected, and we 
loosely refer to it here as the Jeans mass. 

2.2. Models of Primordial Isocurvature Baryon 
Fluctuations 

PIB-likc models may be described by the root- 
mean-square (rms) of the fractional variation in pri- 
mordial baryon density as a function of mass scale 
M, 

(to) (M) = (to) (*„)(£-)-'-*. 

\PbJ rms \Pb ) rms \M N J 

Here Mjv is an abritrary normalization mass scale, 
(<5/3b/ j Ob)(Afjv) is the fluctuation rms on M^r, and n 
is a spectral index. For the primordial abundance 
calculations, however, the knowledge of the statistics 
of the distribution is not sufficient. Rather, we need 
to specify a distribution in baryon density as a func- 
tion of spatial coordinate, pb{x). It is not a unique 
procedure to infer realizations of a distribution in 
baryon-to-entropy from Eq. ||. We follow JF95 and 
consider a Gaussian random variable, A(x). We gen- 
erate stochastic baryon-to-photon distributions, tj(x), 
by the relations 

r l (x)=r !N A 2 (x) , (9) 

or 

r,[x) = mA w {x) . (10) 

Here 77^ is a normalization constant. Note that the 
rj(x) in Eq. ^ and Eq. [l(] are chosen positive definite to 
avoid the introduction of antimatter domains. It has 
been shown in JF95 that the resultant distribution in 
rj is highly non-Gaussian on small scales, in partic- 
ular the transformations in Eq. [)] and Eq. [l^ intro- 
duce phase correlations between the different Fourier 
modes. However, the distribution quickly approaches 
Gaussian character on larger scales. 

We will employ a fluctuation cutoff length scale, 
A c , below which any regions are assumed to be ho- 
mogeneous. As long as we choose A c < e^ookcv how- 
ever, there will be no physical relevance to this scale 
since baryon diffusion prior to the nucleosynthesis 
epoch will provide a natural cutoff scale. Aside from 
any non-Gaussian character of the distribution, the 
statistics of the distribution is then fully described by 
three quantities. These are the effective spectral in- 
dex n in Eq. 0, the rms fluctuation on the cutoff scale, 
(Spb/pb)\ c , and the magnitude of the cutoff scale, A c . 
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2.2.1. Models without Collapse 

We will study three models of small-scale stochas- 
tic baryon inhomogeneity with fluctuations in rj span- 
ning a range from very small scales ~ ^sookeV ^° quite 
large scales well above G?" keV We refer to these mod- 
els as Model 1,2, and 3. They are characterized by 

Model 1 : n « -3, (^j (A c ) = y/2 , (11) 

Model 2 : n pb -3, (A c ) r> 27 , (12) 

Model 3 : n sa -1.5, (A c ) pb 27 . (13) 

In all three models we employ a cutoff-scale A c = 6.8 x 
10~ 5 pc, which is only slightly larger than d^ 00kcV . 
Model 1 is generated by employing Eq. Hand Models 
2 and 3 are generated by employing Eq. |10[ It is seen 
from Eq. || that Models 1 and 2 are chosen to have 
a close to scale invariant spectrum, (Spt/pb) ~ const. 
Note that these models do not incorporate any as- 
pects of structure formation such as the collapse of 
high-density, super-Jeans-mass size regions. The spe- 
cific realizations of the distributions characterized by 
Eq. |ll| - Eq. [l^ which are used to calculate big bang 
nucleosynthesis abundance yields for Models 1-3 are 
shown in the top panels of Figures 1-3. These panels 
show the spatially varying baryon-to-photon ratio di- 
vided by the average baryon-to-photon ratio as a func- 
tion of spatial coordinate. In order to produce these 
distributions we generated 2 x 10 4 Fourier modes. Our 
simulation therefore has a dynamic length scale range 
of 2 x 10 4 . 

It is evident that we examine a one-dimensional 
planar analogue to the full three-dimensional the- 
ory. This is a necessity in order to have sufficient 
length scale coverage to confidently calculate primor- 
dial abundance yields when baryon diffusion is occur- 
ring. Our results are therefore only an approximation 
to a three-dimensional distribution in rj. Even though 
abundance yields in inhomogeneous nucleosynthesis 
have been found to be geometry dependent (Math- 
ews et al. 1990), we expect our one-dimensional cal- 
culations to uncover the trends and magnitudes of 
the effects of baryon diffusion on the nucleosynthesis 
yields. 

There is an additional uncertainty regarding the 
generality of our results for the abundance yields com- 
puted from the specific distributions shown in Fig- 
ures 1-3. Different distributions than those shown in 



Figures 1-3, albeit with the same statistical proper- 
ties, may well give different results for the primordial 
abundance yields. We, however, expect this form of 
"simulation variance" (similar to the variance in simu- 
lations of the formation of large-scale structure) to be 
small in magnitude. Nevertheless, it should be kept in 
mind that particularly for the scale invariant spectra 
Model 1 and 2 there is a small amount of simulation 
variance. 

2.2.2. Models with Collapse 

In the second part of our study we will combine 
Models 1-3 for small-scale stochastic baryon inhomo- 
geneity with large-scale PIB-like fluctuations and the 
collapse of overdense, super-Jeans mass size regions. 
We will refer to these models as Models 1C - 3C. 
We follow JF95 and exclude the abundance yields of 
regions from the abundance average, whenever the 
fractional variations in baryon density exceeds a crit- 
ical value, {Sp/p) > A cr , and when the region's mass 
scale satisfies, M > Mj. As argued in JF95, for 
slightly non-linear, primordial fluctuations, A cr > 1, 
we expect early collapse and, possibly, the formation 
of dark remnants. 

Similarly to the spectrum of small-scale inhomo- 
geneities, and as done in JF95, we will have to in- 
troduce a fluctuation cutoff mass scale for the large- 
scale PIB fluctuations. This is due to the impossibil- 
ity of generating Fourier modes covering the whole 
mass scale range from the very small mass scales, 
M ~ 10 _22 M©, to the mass scales relevant for struc- 
ture formation, M > 3 x 10 5 Af Q . In fact, in or- 
der to calculate the abundance yields for PIB spec- 
tra with small- and large-scale power, we will only 
slightly modify the procedure of JF95. Whereas JF95 
assumed any region with mass M < to be homo- 
geneous, we will assume that any such region's abun- 
dances are well approximated by the abundance yields 
computed from the small-scale inhomogeneous distri- 
bution shown in the top panels of Figure 1-3. 

It is common practice in the field of structure for- 
mation to describe the behavior of fluctuation ampli- 
tudes on length (mass) scale by a primordial power 
spectrum, P(k). In the absence of any significant 
phase correlations the rms fluctuation amplitudes in 
the fractional variations of primordial baryon density 
as a function of mass scale can be derived from the 
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power spectrum by 

/ x \ 2 r k ( M ) 

(J±\ (M)= dink k 3 P(k) . (14) 

V Pb J rms Jo 

In this expression we derive the relationship for k{M) 
from Eq. [j] with A = (27r/fc), corresponding to sharp 
fc-space filtering. The approximate primordial fluctu- 
ation power per unit logarithmic wave vector inter- 
val (e.g. fc 3 / 2 P 1//2 (fc)) which characterizes our models 
of inhomogeneous nucleosynthesis and the collapse of 
overdense regions are shown in Figure 4. Note that 
the figure shows the assumed truly primordial and 
unprocessed power spectra. In particular, any micro- 
physical processing and, for example, different grav- 
itational growth factors for different Fourier modes, 
usually accounted for by a transfer function, are not 
included. 

The behavior of these models on the nucleosynthe- 
sis scale is shown on the left hand (large fc) side of Fig- 
ure 4, whereas the behavior of the models on the large 
scales relevant for collapse is shown on the right hand 
(small k) side of Figure 4. The cutoff for each model 
is shown on the right. In this figure it can be seen 
that the extended mass scale range in our problem 
forced us to resort to rather contrived power spectra. 
In particular, all three models do not have any fluc- 
tuation power on length scales roughly in the range 
Ipc ^ A < IMpc. These power spectra, however, do 
approximate a scenario with a white noise (n = 0) 
large-scale tail, turning over to a scale-invariant spec- 
trum (n — —3) slightly above the Jeans-mass scale, 
and either staying scale-invariant through the scales 
relevant for diffusion during nucleosynthesis (Model 
1C, 2C), or having a second turnover on these scales 
approaching spectral index n = —1.5 (Model 3C). 

Specifically, Models 1C-3C are characterized by the 
following parameters: 

Model 1C : n s r> -3; (A*) = ^2x^2; 

A* = 6.8 x lCT 5 pc , (15) 
»i«0; (jf)(Ai) = V2; 
A^ = 6.8 x lO^/i^Mpc; A cr = 1.5 , 

Model 2C: n s « -3; (A*) = V2 x 27; 

A* = 6.8 x 10" 5 pc , (16) 



n^O; M(X l c ) = 27; 
X l c = 9.5 x 10- 2 /i _1 Mpc; A cr = 2 , 

Model 3C: n s « -1.5; (A*) = y/2 x 27; 

A* = 6.8 x l(T 5 pc , (17) 

m^O; (V)(A<) = 27; 

Aj, = 9.5 x lO^/i^Mpc; A cr = 2 . 

Here the upper lines for each model give the parame- 
ters for the small-scale fluctuations, whereas the lower 
lines give the parameters for the large-scale fluctua- 
tions. Note that the small-scale distribution of Model 
1C (2C, 3C) is identical to the small-scale distribution 
of Model 1 (2, 3). The factors of V2 in the (6p/p)(X s c ) 
have only been written to illustrate that the large- 
scale fluctuation rms has to be added to the small- 
scale fluctuation rms in quadrature. The cutoff scales 
X l c and fluctuation amplitudes {8pbJ Pb)(X l c ) have been 
chosen such that the PIB fluctuation parameters are 
roughly normalized to yield an rms in the fractional 
variation of baryon density on the scale 8ft. _1 Mpc of 
unity at the present epoch (JF95). Here we assumed 
full rcionization. It is evident that the chosen cut- 
off scales for the large-scale PIB-fluctuations satisfy 
X l c > Xj. 

2.3. Numerical Techniques for the Calcula- 
tion of Primordial Abundance Yields 

The small-scale inhomogeneous nucleosynthesis cal- 
culations made use of the nucleosynthesis code devel- 
oped by Kurki-Suonio et al. (1989, 1990a, 1990b, 
1992). The nuclear reaction rates were updated ac- 
cording to Kawano et al. (1993). Regarding weak 
reaction rates, we used the Wagoner (1973) polyno- 
mial for the weak reaction rates, with r„ = 887.0s. 
We included all recently determined corrections to 
this rate as relevant (Dicus et al. 1982; Kernan & 
Krauss 1994; Seckel 1993; Kernan 1993). We also cor- 
rected for effects of finite spatial resolution and time 
step effects. The net correction factor to Y p from all 
of these effects was only -0.0012. (Note that the cor- 
rection + .0017 found by Kernan (1993) by using a 
shorter time-step does not apply to our code since the 
code used here is independent from the one used by 
Kernan.) We use all the relevant diffusion coefficients 
for protons and neutrons as given by Jedamzik & 
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Fuller (1994). However, we do not incorporate the 
effects of late-time hydrodynamic expansion and neu- 
trino cooling of inhomogeneities into our calculations 
(Jedamzik, Fuller, & Mathews 1994) which renders 
our calculated ( 7 Li/H) ratios somewhat higher than 
when hydrodynamic expansion is included. 

A key challenge in the present work was the neces- 
sity to utilize a very large spatial grid (40 000 zones). 
It was thus essential to optimize the speed of compu- 
tations. We therefore included only the 8 most im- 
portant light isotopes in our reaction network. The 
reactions leading to isotopes heavier than A= 7 were 
included as sinks in the reaction flow. Since there is 
little flow from heavier isotopes back to lighter, this 
is an excellent approximation for the present applica- 
tion. The most time consuming part of the nucleosyn- 
thesis calculation is the inversion of an 8 x 8 reaction 
matrix which must be performed at each zone at each 
time-step. Since this inversion does not easily vector- 
ize, we found it best to run on a fast scalar machine. 

3. Results 

In this section we present the results of our inhomo- 
geneous primordial nucleosynthesis calculations. Fig- 
ures 1-3 show the employed random distributions in 
r](x)/(rj) for the small-scale baryon inhomogeneities 
of Models 1-3 (and Models 1C-3C) in the top pan- 
els. Nucleosynthesis yields for the various isotopes 
and models are given along the side panels of Figures 
1-3. The three panels on the left-hand-side of the fig- 
ures give the 4 He mass fraction, Y p , ( 2 H/ 1 H) -number 
ratio, and ( 7 Li/H) -number ratio as a function of av- 
erage baryon-to-photon ratio r\ for the models with- 
out collapse. The three panels on the right-hand-side 
of the figures give the abundance yields for the same 
light elements as a function of baryon-to-photon ratio 
for the models with collapse. For reference the frac- 
tional contribution of baryons to the critical density, 
fib , can be obtained from the baryon-to-photon ratio 
r] by fl b = 0.0112(77/3 x 10- w )h- 2 . 

3.1. Models without collapse 

For the panels on the left-hand-side the three lines 
correspond to: Solid lines-the results of a homoge- 
neous big bang nucleosynthesis (hereafter; BBN) cal- 
culation at the same average 77; Dotted lines- the 
results of an inhomogencous BBN calculation of the 
distribution shown in the top panel of the figure but 
without taking into account of the effects of baryon 



diffusion; and Dashed lines- the results of our inho- 
mogencous BBN calculations which include the effects 
of baryon diffusion. 

In the left-hand-side panels of Figure 1 we show 
nucleosynthesis yields for Model 1 . It is seen that the 
abundance yields of the inhomogeneous BBN Model 

1 are not very different from the yields of a homo- 
geneous BBN scenario at the same average 77 (with 
the possible exception of 7 Li). This is mainly be- 
cause the fluctuation amplitudes in Model 1 are rel- 
atively small. The trends in Model 1 can be best 
understood when compared to the results of Model 

2 which are shown in the left-hand-panels of Figure 
2. Note that both Model 1 and 2 have an approx- 
imately scale-invariant spectrum (n w —3) but that 
Model 2 has much larger fluctuation amplitudes than 
Model 1. It is evident, particularly from the results 
of Model 2, that for close to scale-invariant spectra 
baryon diffusion does not significantly change the cal- 
culated abundance yields. This can be inferred by 
comparing the calculated abundance yields of the in- 
homogeneous distribution where baryon diffusion is 
accounted for (dashed line) to those where baryon 
diffusion is neglected (dotted line). This trend is also 
observed in Model 1 (Figure 1 ) . That baryon diffusion 
is comparatively unimportant can be understood by 
inspection of the baryon distributions in Figure 1 and 
2. The bulk of the baryon number in both models is 
included in fairly extended high-density regions with 
only relatively small variations in 77 around the ap- 
proximate average baryon-to-photon ratio of these ex- 
tended regions. It is known that the effects of baryon 
diffusion are maximized when significant fluctuations 
exist on scales comparable to the neutron diffusion 
scale during the epoch of primordial nucleosynthesis 
(Applegate, Hogan & Scherrer 1987). The effects of 
baryon diffusion on Models 1 and 2 arc small since the 
typical length scale of extended high-density regions 
is much larger than the neutron diffusion length scale. 

For most 77 values the 7 Li yields are substantially 
higher for the inhomogeneous distribution than the 
corresponding Li yields in a homogeneous BBN sce- 
nario. This is a generic feature of inhomogeneous pri- 
mordial nucleosynthesis (Alcock et al. 1987). Be- 
cause of this, inhomogencous BBN scenarios are only 
compatible with observations if there is moderate to 
significant depletion of 7 Li in low-metallicity Pop II 
halo stars. In fact, there may be evidence for 7 Li 
depletion in Pop II halo stars (Pinsonncault ct al. 
1992). However, the situation remains unresolved 
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(Schramm & Mathews 1996). 

If one allows for the possibility of 7 Li depletion, and 
assumes that the primordial ( 2 H/ 1 H) -ratio is as large 
as may have been inferred from some (but not all) 
observations of deuterium in Lyman-limit clouds at 
high redshift (Songailia et al. 1994; Rugers & Hogan 
1996; but sec also Tytler, Fan, & Buries 1996; Buries 
& Tytler 1996), Model 1 is marginally consistent with 
observations for some range in r/. In contrast, Model 
2 has no range in r\ consistent with the observation- 
ally inferred light-element abundances. We conclude 
that the fluctuation amplitudes of stochastic, approx- 
imately scale-invariant baryon inhomogeneities, are 
severely constrained by the possible overproduction 
of 4 He and/or 7 Li during the epoch of big bang nu- 
cleosynthesis even when the effects of baryon diffusion 
are included. 

In Model 3 we assumed a fluctuation spectrum 
with spectral index n w —1.5 and with significant 
fluctuation amplitudes. In contrast to Models 1 and 
2 this model is characterized by decreasing fluctua- 
tion amplitudes with increasing length scale as evi- 
dent from the r\ distribution shown in the top panel 
of Figure 3. In a similar fashion to Figures 1 and 2 the 
left-hand-side panels of Figure 3 show the abundance 
yields of Model 3. 

The results of Model 3 are, perhaps, somewhat 
startling. It is evident from comparison of the dot- 
ted and dashed lines in Figure 4 that for the param- 
eters employed in Model 3 baryon diffusion signifi- 
cantly changes the abundance yields. In particular, 
the 4 He mass fraction in Model 3 is very close to that 
of a homogeneous BBN scenario and much smaller 
than the calculated 4 He mass fraction for the inho- 
mogeneous distribution when baryon diffusion is ne- 
glected. For most of the range in rj 7 Li yields are 
substantially smaller and deuterium yields are larger 
than those in inhomogencous BBN scenarios without 
diffusion. 

The effects of baryon diffusion are important since 
significant fluctuation amplitudes exist on the scale of 
the neutron diffusion length during the epoch of pri- 
mordial nucleosynthesis. Fluctuations on larger spa- 
tial scales, which are much less affected by neutron 
diffusion and which may grossly overproduce 4 He and 
7 Li, are substantially smaller than those on the neu- 
tron diffusion length due to the large spectral index 
of the distribution. In other words, the average rj in 
extended regions is close to the average baryon-to- 
photon ratio for the whole distribution. 



Assuming high 2 H and significant 7 Li depletion in 
Pop II halo stars there is some range for rj in Model 3 
which can agree with all the observationally inferred 
primordial abundance yields. This is quite surprising 
since it implies that random variations in the baryon- 
to-entropy with rms amplitude as large as (dpb/pb) ~ 
30 on c^pokev not ruled out by the observations 
of 2 H, 4 He, and 7 Li in low-metallicity environments. 
Nevertheless, if in future it is shown convincingly 
that either the primordial ( 2 H/ 1 H) P w 3 x 1CP 5 is 
low or the primordial ( 7 Li/H) p w 2 x 1CP 10 is low, 
then significant baryon-to-entropy fluctuations during 
the epoch of primordial nucleosynthesis are probably 
ruled out. This applies for both, inhomogeneities in 
a nearly regular lattice (such as might be generated 
by a first-order QCD phase transition), or stochastic 
inhomogeneities of arbitrary statistical properties as 
may be generated during inflation. 

3.2. Models with collapse 

We now turn our attention to the results of Models 
1C, 2C, and 3C which incorporate the collapse of over- 
dense, large-scale regions. The abundance yields from 
such overdense, super- Jeans-mass size regions will be 
excluded from the abundance average since these re- 
gions either may form black holes, brown dwarfs, 
and/or stellar remnants, or might undergo sufficient 
early stellar processing and chemical evolution which 
would render the regions unsuitable for the determi- 
nation of primordial abundances (JF95). The results 
for the primordial abundances of Model 1C, 2C, and 
3C arc presented in the right-hand-side columns of 
Figure 1, 2, and 3, respectively. Each panel in these 
figures shows the results for: Dotted lines - a model 
which does not have any sub-Jeans-mass size fluctua- 
tion power as in JF95; Dashed lines - a model which 
includes fluctuations on scales below the Jeans length 
and explicitly treats diffusion as specified in Section 
2; and for comparison, Solid lines - the results of a 
standard homogeneous BBN calculation. 

As already discussed in detail in JF95, models 
which have a fluctuation cutoff above the Jeans mass, 
and which have high efficiencies for collapse of over- 
dense regions, generally yield substantially lower 4 He 
mass fractions Y p , higher ( 2 H/ 1 H) -number ratios, 
and comparable to higher ( 7 Li/H) -number ratios 
than a homogeneous BBN scenario. The modifica- 
tions to these results by the introduction of small- 
scale power on scales below the Jeans scale studied 
here are evident from the dashed lines. In general, 
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small-scale fluctuations enhance 4 He and 7 Li yields 
over those yields from PIB models without small-scale 
power. This is easily understood since regions fluctu- 
ating in -q mostly yield larger Y p and ( 7 Li/H) ratios 
(even when baryon diffusion is included) than regions 
which are homogeneous in 77 (as apparent from the 
left-hand-side panels of Figures 1-3). 

If PIB models with collapse and small-scale fluctu- 
ations even exceed the 4 He and 7 Li yields of a stan- 
dard homogeneous BBN scenario largely depends on 
rms fluctuation amplitude. Whereas Model 1C is 
marginally compatible with observationally inferred 
primordial abundances Models 2C and 3C, which 
have large rms amplitudes in 77, either overproduce 
7 Li and/or 4 He for large 77 or overproduce deuterium 
for small 77. It is evident that none of the models 
which include small-scale fluctuations can agree with 
observations for large average baryon-to-photon ra- 
tios (77 ^ 3 x 10 -10 ) mainly due to an overproduction 
of 7 Li. This is not the case for PIB models which are 
homogeneous on scales below the Jeans scale (JF95). 
This implies that the introduction of small-scale fluc- 
tuations into stochastic PIB models considerably nar- 
rows, or even excludes, the ranges of 77 which may 
have agreed with observationally inferred primordial 
abundances for models without small-scale fluctua- 
tions. 

4. Conclusions 

We have investigated the primordial nucleosynthe- 
sis process for models which include random varia- 
tions in the baryon-to-photon ratio. Our study ex- 
plicitly accounts for the effects of baryon diffusion. 
Analyzing three different stochastic distributions in 77, 
which have different fluctuation amplitudes and spec- 
tral indices, we have found that the statistical prop- 
erties of stochastic baryon-to-photon inhomogeneities 
are severely constrained by the possible overproduc- 
tion of either 7 Li and/or 4 He for large average 77 or 
deuterium for small average 77. In particular, we have 
found that the rms amplitudes of scale-invariant 77 
distributions, characterized by almost equal fluctua- 
tion amplitudes on all scales, are constrained to be 
not larger than a few. 

In contrast, the constraints on the rms amplitudes 
of distributions characterized by decreasing fluctua- 
tion amplitudes with increasing mass scale are much 
weaker. In such models rms amplitudes as large as 
~ 30 on the neutron diffusion scale at weak freeze-out 



are currently not ruled out by observationally inferred 
primordial abundance limits. Nevertheless, these con- 
straints may be considerably tightened if it is ever 
convincingly shown in the future that either depletion 
of primordial 7 Li in low-metallicity halo stars is neg- 
ligible or that the primordial deuterium-to-hydrogen 
ratio is low ( 2 H/ X H) - 3 x 1CT 5 . 

In the second part of our study we have considered 
stochastic baryon inhomogeneity extending from the 
very small baryon diffusion scale during the epoch of 
primordial nucleosynthesis to very large scales above 
the post-recombination Jeans length, and the early 
collapse of overdense super- Jeans-mass size regions. 
In agreement with JF95 we have found that signifi- 
cant fluctuation amplitudes on scales below the Jeans 
scale considerably narrow (or exclude) the parameter 
space of such PIB-like models which may agree with 
observationally inferred primordial abundance limits 
otherwise. We conclude that PIB fluctuations may 
only be consistent with observationally inferred pri- 
mordial abundance limits for large fib if either fluc- 
tuations are highly phase correlated (Gnedin et. al. 
1995) or if there exists a fluctuation cutoff on scales 
above the Jeans scale (JF95). 

This work was performed in part under the aus- 
pices of the US Department of Energy by the Lawrence 
Livermore National Laboratory under contract num- 
ber W-7405-ENG-48. Work at University of Notre 
Dame was supported in part by DOE Nuclear The- 
ory grant DE-FG02-95ER40934. Computation was 
done in part at the Center for Scientific Computing 
(Finland) . 
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Fig. 1. — Baryon-to- photon distribution(top panel) 
and nucleosynthesis yields for Model 1 (left) and 
Model 1C (right side panels). A value of unity on the 
ordinate of the top panel corresponds to 2 x 10 4 A C for 
Model 1 (or 2 x 10 4 A^ for Model 1C), where A c (A*) is 
specified in the text. The side panels give the calcu- 
lated 4 He mass fraction Y p , ( 2 H/ 1 H) -number ratio, 
and ( 7 Li/H) -number ratio as a function of the aver- 
age r\. The left-hand side panels give results for the 
model without collapse (Model 1). The dashed lines 
are for the inhomogeneous distribution when baryon 
diffusion is included. The dotted lines are for the 
inhomogeneous distribution when baryon diffusion is 
neglected. For comparison, the solid lines are abun- 
dance yields for a homogeneous BBN calculation. The 
panels on the right-hand-side show the abundance 
yields for models with collapse of high density regions 
(Model 1C). The dashed lines are for an inhomoge- 
neous distribution with collapse when small-scale fluc- 
tuation as shown in the top panel are present. The 
dotted lines are for an inhomogeneous distribution 
with collapse when there are no small-scale fluctu- 
ations. The solid lines are for a homogeneous BBN 
calculation. Note that the large-scale fluctuation part 
of Model 1C is not shown in this figure. 



Fig. 2. — Same as Figure 1, but for Model 2 and 
Model 2C. 



Fig. 3. — Same as Figure 1, but for Model 3 and 
Model 3C. 



Fig. 4. — Illustration of the primordial fluctua- 
tion power per unit logarithmic wave vector inter- 
val P2(k)kl as a function of wave vector k for the 
models employed in the nucleosynthesis calculations. 
The steeply sloping lines on the far left are a sim- 
ple extrapolation of a PIB model with n = 0. This 
extrapolation is disrupted as indicated for the three 
models. At the nucleosynthesis scale, three possible 
fluctuation spectra are considered as described in the 
text. 
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